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I. INTRODUCTION 



Development of ab initio theories for neutral open- 
shell atoms is a difficult task since the interaction be- 
tween electrons of an open shell is strong and cannot be 
treated perturbatively. Nevertheless, some progress in 
two- valence and even three- valence electron atoms has 
been achieved with the combination of configuration- 
interaction (CI) method and many-body perturbation 
theory (MBPT) @, §, {§. Particle-hole states of closed- 
shell atoms have additional difficulty that conventional 
perturbation theory does not converge for hole states. We 
solved the convergence problem by modifying denomina- 
tors H [|. As a result, we were able to achieve agreement 
with experiment for neon energies of J=l odd excited 
states and oscillator strengths (averaged over many mea- 
surements) of transitions to the ground state. In this pa- 
per, we extend our application of the CI+MBPT method 
to neon transition between excited states. If we succeed, 
our understanding of neon atom and low-Z neonlike ions 
will be substantially improved 

Fairly accurate (about 5%) measurements of many 
transition rates between excited states are available, pro- 
viding important tests of theory. In addition, semiem- 
pirical calculations can be compared with our calcula- 
tions. For example, many transition rates along the 
neon isoelectronic sequence were calculated by Hibbert 
et al. |J with a general configuration-interaction code 
(CIV3) In calculations a few parameters were ad- 
justed to fit experimental energies. However, even after 
such adjustments, the results still disagreed significantly 
with other semiempirical calculations by Seaton and 
experiments. The latter theory was more successful, giv- 
ing results in close agreement with experiments. No pure 
ab initio theory, as far as we know, was (successfully) 
applied previously to calculations of transitions between 
neon excited states. For transitions to the ground state, 
elaborate ab initio calculations exist (Avgoustoglou and 
Beck M), but agreement with experiment for an oscilla- 
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factory. These calculations are more effective in heavier 
noble-gas atoms, where the agreement with experiment is 
achieved for Ar, Kr, and Xe atoms. In neon, on the other 
hand, our new calculations with CI+MBPT method, give 
results close to the average experimental values. 

It is a well-known fact that neon transitions between 
excited states are sensitive to the accuracy of fine- 
structure splittings. Semiempirical theories avoid this 
difficulty by introducing and adjusting several parame- 
ters to match energies of multiplets as precisely as possi- 
ble. For example, using quantum defect method, Seaton 
|Q was able to obtain very small root-mean square devia- 
tions for energies. As a result, he also was able to obtain 
transition oscillator strengths that agree well with ex- 
periment. Getting accurate fine structure intervals with- 
out parametric adjustments is a challenging task. We 
will demonstrate in this paper that CI calculations cor- 
rected with second-order MBPT give energies and os- 
cillator strengths for transitions between excited states 
with the precision comparable to the precision of best 
semiempirical calculations. 

The transition data in neon and other noble gases are 
needed for plasma physics and studies of discharges with 
many industrial applications in lamps and gas lasers. The 
opacity project |J is another motivation behind many 
calculations (one example is given in Ref. fl0|| ) in neon- 
like ions. Understanding of neon atom can be beneficial 
for the development of atomic structure methods which 
are needed for many applications. One important ap- 
plication of atomic structure is calculations of parity- 
nonconservation amplitudes in heavy atoms with one or 
a few valence electrons which require a clear understand- 
ing of correlation effects in these atoms. One-valence- 
electron MBPT has similar convergence problems to the 
hole MBPT after a core is excited. Modification of de- 
nominators according to our prescription might be one 
key to the solution of a puzzling problem that third-order 
energy in Cs agrees worse with experiment than second- 
order energy. Calculations of electron dipole moments 
in particle-hole atoms is another, more direct application 
of our particle-hole theory. Furthermore, the CI+MBPT 
method and convergent hole perturbation theory can be 
generalized for more complicated atoms with more than 
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one particle or hole and properties of these atoms can be 
explored beyond Hartree-Fock approximation. 

In this paper, first we will briefly describe our method 
of calculations (more details are given in Ref. [|| ||); then, 
we will compare CI+MBPT and experimental energies 
for J=l and J=2 even states. This comparison gives an 
estimate on accuracy of our wave functions. Next, we will 
show our results for transition line strengths. Finally, we 
compare our theory with other semiempirical calculations 
and experiments. 



II. CI+MBPT CALCULATIONS 

A. Energies and oscillator strengths for J=l odd 
neon states 



TABLE I: Calculations of neon energy levels for J=l even 
states. In the third and forth columns "CI-8" and "CI-32" 
mean that in our calculations the size of the CI matrices were 
8 and 32, respectively. A is the difference between theoretical 
and experimental energies. All energies are in atomic units. 



J=l even NIST CI-8 CI-32 A shifted A 

P3/ 2 3p 0.67551 0.6687 0.6690 0.0065 0.0005 

Pg/ 2 3p 0.68400 0.6789 0.6787 0.0053 -0.0007 

Pi/ 2 3p 0.68696 0.6820 0.6817 0.0052 -0.0008 

Py 2 ?>p 0.68818 0.6833 0.6830 0.0052 -0.0008 

p-^p 0.74048 0.7337 0.7339 0.0066 0.0006 

p-^p 0.74274 0.7370 0.7367 0.0060 0.0000 

p-^p 0.74567 0.7398 0.7396 0.0061 0.0001 

Pw 2 4p 0.74590 0.7402 0.7399 0.0060 0.0000 



The Rayleigh-Schrodinger variant of second-order 
MBPT, given in Ref. (T^|, has low accuracy for neon 
(does not improve lowest order approximation), there- 
fore we developed a fast convergent variant of MBPT 
j| |[. This perturbation theory can be understood 
from couple-cluster single-double equations Sim- 
ply put, we modify some denominators in the pertur- 
bation terms to take into account the strong interac- 
tion between a hole and a core electron or between core 
electrons non-perurbatively. The advantage of this ap- 
proach compared to the couple-cluster method of Q 
is simplicity and speed of calculations. With our fast 
convergent MBPT method, we are able to improve the 
accuracy of hole energies and fine-structure splittings 
of light neonlike ions already after adding second-order 
MBPT corrections. Apart from Coulomb correlation 
corrections, the Breit magnetic interaction is also in- 
cluded, but small frequency-dependent Breit, quantum- 
electrodynamic, reduced-mass, and mass-polarization 
corrections are omitted (the analysis of these small cor- 



rections is given in Ref. 12 



To calculate particle-hole energies, we construct a 
model CI space [pT) , compute effective Hamiltonian in 
this space, which also includes second-order MBPT cor- 
rections, and solve an eigen value problem. Along with 
energies we obtain wave functions, which were used to 
calculate oscillator strengths for transitions to the ground 
states H |j. The energies of neon particle-hole J=l odd 
states and oscillator strengths were in very good agree- 
ment with experiment after using a relatively small CI 
space, 52. Pure ab initio energies differed from experi- 
mental energies by 0.0069 a.u., but after subtraction of 
the systematic shift (which does not make much differ- 
ence in transition calculations), the agreement was im- 
proved to the level of 0.0001 a.u. for almost all states. 
We will use the same wave functions for J=l odd states in 
our calculations of transitions from J=l odd to J=l and 
J=2 even states. The accuracy of energies of even states 
involved in the transitions will be illustrated below. 



TABLE II: Calculations of neon energy levels for J=2 even 
states. The size of the CI matrix is 32. A is the difference 
between theoretical and experimental energies. All energies 
in atomic units 



J =2 even 

P3/2 3 P 
P3/2 3 P 
Pl"/2 3 P 
P3/2 4 P 
P3/2 4 P 
Pl~/2 4 P 



NIST CI-32 



A shifted A 



0.68265 0.6775 0.0051 -0.0007 

0.68489 0.6800 0.00487 -0.0009 

0.68736 0.6825 0.00485 -0.0010 

0.74222 0.7363 0.00588 0.0001 

0.74285 0.7372 0.00565 -0.0001 

0.74591 0.7401 0.00578 -0.0000 



B. Energies for J=2 and J = l even neon states 

In this section we present our calculations of energies 
for even states. We use the same formalism and numeri- 
cal method as in Ref. |||. The spline cavity is chosen 
80 a.u., the number of splines is chosen 40, and the max- 
imum orbital momentum is chosen 5. For excited states, 
a V^" 1 Hartree-Fock (HF) potential basis (see [||) is 
built from HF V N spline basis by diagonalization of a 
one-electron Hamiltonian to take into account the major 
part of the interaction of an excited electron with a hole. 
Such procedure speeds up the convergence of CI and re- 
duces uncertainty in denominators of a perturbation the- 
ory. A V N ~ X HF basis used in Ref. jjJI was constructed 
by solving differential HF equations and gives similar re- 
sults, but our basis does not require rewriting HF code 
and therefore is more convenient. The calculations of 
energies are shown in two tables. 

In Table || we compare with experiment our theoreti- 
cal energies of 3=1 even states. Energies calculated in 
a model space of size 32 agree better with experiment 
than energies calculated in smaller size (8) model space. 
CI space 32 can be considered as optimal since a larger 
number of configurations does not improve much accu- 
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TABLE III: Transitions between neon excited states, from 
J=l odd to J=l even. In notations, p"^ and Py 2 are no ^ e 
states, and particle states are standing immediately to the 
right. For unique specification, experimental wavelengths are 
also provided. The sizes of model space for even states are 
given, but for J=l odd states the size is 50, the same in all 
cases. The experimental NIST and theoretical line strengths 
(in columns "CI-8" and "CI-32") are expressed in a.u.; A 
denotes the relative deviations of the theoretical line strengths 
from the experimental line strengths 

Transition A, A NIST Ac. CI-8 CI-32 +RPA A 

P3/2 3s -P3/2 3 P 
p-^s-py^p 

Pi/2 3s -Pl/2 3 P 

~/2 4 P 

J 3/2' 



-P3/2' 

-P3/2 3 P 
-P3/2 3 P 



^3/2 
PS/2 
Pl/2 
PV2 3 S-P 3 /2 4 P 
Pl>2 3S -Pl/2 3 P 
Pl/2 3s -P372 3 P 
Pl/2 3s -Pl>2 3 P 
Pl/2 3S -Pl/2 4 P 
Pl/2 3S -Pl/2 4 P 



6385 


12.4 


B- 


12.5 


12.6 


12.2 


2% 


6032 


1.82 


B- 


1.83 


1.88 


1.81 


1% 


6130 


0.23 


B- 


0.32 


0.23 


0.22 


4% 


3502 


0.076 


D 


0.114 


0.102 


0.084 


11% 


7247 


5.27 


B- 


5.99 


5.93 


5.70 


8% 


8085 


0.094 


B 


0.121 


0.107 


0.103 


10% 


3687 


0.029 


D 


0.038 


0.040 


0.032 


10% 


6601 


9.88 


B- 


9.16 


9.25 


8.92 


10% 


7026 


0.971 


B 


0.825 


1.01 


0.982 


1% 


6719 


9.75 


B- 


11.0 


10.9 


10.5 


8% 


3595 


0.045 


D 


0.070 


0.065 


0.053 


18% 


3601 


0.030 


D 


0.050 


0.045 


0.036 


21% 



racy. Note that the size of optimal CI space depends on 
the choice of a starting potential since inadequate ini- 
tial approximation is corrected by diagonalization of the 
CI matrix. An 1% systematic shift is present, which can 
be attributed to inaccuracy of hole energies; however this 
shift does not affect much accuracy of transition rates and 
is relatively unimportant. More important for weaker 
transitions is the fact that after subtracting this shift, 
we obtain very small residual deviations, which bode well 
for the accuracy of singlet-triplet mixing coefficients and 
transition amplitudes. 

Similar agreement of energies is obtained for 3=2 even 
states in Table [u|. Again, we have almost the same sys- 
tematic shift, and after its subtraction, only a very small 
difference between experimental and theoretical energies 
remains. Therefore, we have reason to expect good preci- 
sion for transitions between the excited states which are 
considered next. 



TABLE IV: Line strengths (a.u.) for transitions between neon 
excited states from J=l odd to J=2 even. Abbreviations are 
the same as in Table 
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P3/2 
P3/2 



3s- 
3s- 



Transitions 

=l 3 ~ PV2 3 P 
P 3 /2 3 P 
P3/2 3 P 
P3/2 4 P 
Pl/2 3S -P3/2 4 P 

Pi72 3s -Pl/2 3 P 
•*• ~~2 3 P 



A, A 



NIST Acc. CI-32 +RPA A 



p 3/ , 2 3s 
" 3s- 



P3/2 



Pl/2 



3s- 

Pl/2 3s - 

Pi72 3s - 



P3/2' 
Pl/2 4 P 
P. ; 



^3/2 3 P 



6097.8507 


10.1 


C+ 


10.2 


9.76 


3% 


6508.3259 


20.4 


B- 


21.1 


20.2 


1% 


6306.5325 


2.57 


B- 


2.54 


2.62 


2% 


3516.1960 


0.074 


D 


0.106 


0.077 


4% 


3702.2783 


0.028 


D 


0.034 


0.029 


4% 


6680.1202 


17.1 


C+ 


16.9 


17.1 


0% 


6931.3788 


14.3 


B- 


16.6 


15.3 


7% 


3594.5516 


0.11 


D 


0.14 


0.11 


0% 


7175.9155 


2.62 


B- 


2.51 


2.61 


0% 



TABLE V: Comparison of experimental and theoretical os- 
cillator strengths of neon. We adopt compact notations of 
Ref. j?|: p t q t ri -p f q f r f , where p = 2xj a ,q = 2xK,r = J, 
i stands for initial, and / stands for final. The hole angu- 
lar momentum j a is coupled with the orbital momentum of 
the excited electron I to give the angular momentum K in 
the intermediate coupling scheme; K is coupled with a spin 
of the excited electron to give a total angular momentum of 
particle-hole state J. For compete specification, nl quantum 
numbers are also provided. Brackets denotes powers of 10. 



Transitions A, A NIST 
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Ref.]il| 



Current Ref. 
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Previously, we calculated oscillator strengths for the 
transitions to the ground state; the formula for excited 
state transitions is different. The final expression for the 
coupled reduced matrix element after angular reduction 
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in j-j relativistic basis has the following form: 



(F\\Zj\\I) = y/[J F ] [J I ]C F (a'v')C I (a,v)x 

{ _ 1)J+ J I+ja+jvl | h J f\ Sa/a{v '\\Zj\\v) + 
{ Jv' Ja Jv ) 

(-l) Jrt " +3 " +1 (j f i F )s v , v (a\\Zj\\a')\ 

^ Jv Ja' Ja J J 

We use standard notations of relativistic MBPT meth- 
ods, see for example fl4|| . Configuration weights of the 
final Cf{o!v') and the initial Ci(a,v) states are obtained 
in CI calculations. Note that this formula is similar to 
but different from that for the transitions between two- 
particle states |l4|] even if we neglect small hole-hole ma- 
trix elements. However, it is possible to modify wave 
functions to use two-particle matrix elements, which can 
be convenient if the program for calculations of two- 
particle matrix elements are available. 

Tables [II and IV show our results of calculations for 
J=l-J=l and J=l-J=2 excited state transitions. Many 
precisely measured (5% level) neon transition rates pro- 
vide an imp ortant test of accuracy of our calculations. 
In Table Ell the calculations are done with two numbers 
of configurations to show improvement in precision for 
the larger number of configurations. Important random- 
phase approximation (RPA) corrections are included by 
replacing "bare" matrix element with "dressed" , RPA 
matrix elements (the replacement for two-valence elec- 
tron atoms was implemented in [§, ||). Including RPA 
corrections needs some care because a hole state present 
when the transition between excited states occurs leads 
to the convergence problem and low accuracy of the reg- 
ular RPA corrections. Our standard cure is to mod- 
ify denominators by subtracting radial Slatter integral 
R^{abab) thus taking into account the monopole inter- 
action of a hole with a core electron. This modification 
of denominators approximately halves RPA corrections. 
Adding divided by 2 normal RPA corrections, we esti- 
mated the level of these corrections and found that they 
are important and improve agreement with experiment. 
In neon, the corrections constitute a few percent part of 
a total matrix element, but in heavier noble gas atoms, 
they are even larger and more important. In the last col- 
umn, we place our best values calculated in CI-32 model 
space with appropriate modified-denominator RPA cor- 
rections. 

The deviation from experiment is consistent with the 
experimental accuracy (for example, the accuracy of class 



"B" is in the range 5% and the deviation from theory is of 
the same magnitude). For transitions which have exper- 
imental accuracy of classes "B-" and "C+" , the theory 
is as accurate as or even more accurate than experiment, 
but for the class "D" , the theory is definitely more ac- 
curate. Still the accuracy of theory is not the same for 
all transitions, since some suppressed transitions owing 
to cancellation could be more sensitive to fine-structure 
splittings. 

Our final table (Table |v|) contains our best values of 
oscillator strengths with the RPA corrections for com- 
parison with experiments and other theories. Oscillator 
strengths are calculated from line strengths S and tran- 
sition energies u> in atomic units. 
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Due to a large number of measurements and calculations 
we restricted ourselves to comparison with results from a 
few sources, which contain further references; for exam- 
ple, references to many experiments and comparison with 
several measurements are given by Bridges and Wiese (ll| 
where the authors also estimated uncertainties of their 
experiment to be about 7% and of the others shown in 
their comparison table to be in the range 10-50%. The 
experiments seem to be in good agreement with each 
other. NIST data derived from various sources are in 
close agreement with values given by Bridges and Wiese 
p3[ . The agreement with experiment of the semiempir- 
ical theory by Seaton is similar to the agreement of 
our ab initio theory. Calculations performed by Hibbert 
et al. |^| agree worse; for example, for the 6130 A transi- 
tion, an experimental value is 0.0114 or 0.0119, our value 
is 0.0109, but the value in Ref. [g is 0.0249. Overall 
agreement of theories and experiments is quite normal. 



III. CONCLUSIONS 

In this paper, we have applied CI+MBPT theory for 
particle-hole states of closed-shell atoms to calculations 
of transitions between excited states of neon. A difficulty 
that the hole energy has poor convergence is overcome 
with modifications of denominators in MBPT. Good pre- 
cision for particle-hole states is illustrated for many en- 
ergy levels of neon. Apart from energies, our theory is 
tested in calculations of line strengths. Agreement with 
experimental values is achieved. 
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